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1. INTRODUCTION 

The operation of electrical power systems (EPS) is commonly assessed through load flow, voltage 
stability, and short-circuit studies, among others. This is done with the help of software programmed with 
electromagnetic equations and mathematic models that simulate the behavior of the main components (i.e., 
generators, load, line, transformers, and machines) [1], [2]. Updating the load models commonly used in 
simulation programs of EPS is an imperative need today due to the increase in emerging smart grid technologies 
such as distributed generators, electric vehicles, and static power converters, among others. Modeling this type of 
load is difficult due to its time-varying and weather-dependent nature, as well as the lack of measurements [3]. 

The three-phase asynchronous motors (TAM) are one of the most important loads in EPS since it is 
estimated that there are around 300 million of this equipment installed in industries [4], representing the 
consumption of around 43 to 46% of energy electricity worldwide and between 60 and 70% of energy in 
industries [5]-[7]. Due to the importance of TAMs, they have been the subject of numerous studies focused 
mainly on the improvement of the technology to increase the level of efficiency [8], [9], on control strategies 
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in the TAMs drive system to guarantee maximum operating efficiency [10]-[12] and on methods for 
estimating TAM operational efficiency [13], [14] 

The importance of improving the modeling of TAMs in the simulation programs of EPS was 
evidenced by the blackout in the United States that occurred in 1996 and 2000 when millions of North 
Americans were without electricity for several hours [15]. In a study led by the "western electricity coordinating 
council," it was shown that the system in the affected area did not include the modeling of various TAMs and 
electronic loads, therefore, the behavior of the EPS could not be predicted correctly [16]-[18]. 

Normally, the TAM is represented as a constant power load, or constant impedance [3], [19]. 
Assuming constant active power may be acceptable in some cases but ignoring the variation of reactive 
power with voltage in TAM leads to significant errors. Neither has the behavior of the power consumed by 
the TAM been analyzed in the event of frequency deviations from its nominal value. Normally, in frequency 
stability studies, concentrated models of the loads are used, for which it is necessary to determine the so- 
called load-damping constant [2]. There are no concrete proposals in the literature to estimate this parameter 
in predominantly TAM loads. 

According to CIGRE [20], there are two types of load models depending on the way to obtain their 
parameters. One way is based on measurements made in situ and, another way is based on the knowledge of 
the components of the loads and the development of models estimating the corresponding parameters. In the 
case of the TAM, the method based on measurements is impracticable because it requires laboratories where 
a considerable number of TAMs can be tested, both low and medium voltage and of different power and 
number of poles. 

This paper aims to propose a new model of the TAM that allows its use in the analysis of EPS. The 
proposed model allows the analysis of the influence of the variation of voltage and frequency deviation, as 
well as the type of the mechanical load driven by the motor, on the active and reactive power consumption of 
the motors. The inclusion of the consideration of the type of mechanical load driven by the motor in the 
composite model is a novel aspect that can contribute to the improvement of the EPSs simulation software. 
Section 2 describes the state of the art of developed models of TAMs as load in the EPS. Section 3 explains 
the model with the equations and the sequential procedure for its application, while the explanation of the 
model is presented in section 4. In section 5 the conclusions of the article are discussed. 


2. RESEARCH METHOD 
2.1. State of the art of developed models of TAMs as load in the EPS. 

Balanathan et al. [21] the approximate equivalent circuit shown in Figure 1 is proposed. In this case, 
the stator resistance and magnetization reactance are not considered. Neglecting the stator resistance does not 
generate a significant error, but not considering the magnetization reactance implies simplifying the main 
consumer of reactive power, which is unacceptable. 

Where: P is the active power, Q is the reactive power, I is the current, V; is the voltage per phase, X; is 
the stator leakage reactance, X2 is the rotor leakage reactance, R2 is the rotor resistance, E is the electromotive 
force induced in the rotor, and s is the slip. Choi et al. [22] it is proposed to incorporate the equivalent circuit, as 
a model, in power system studies, as it appears in Figure 2, giving rise to the so-called "ZIP model with TAM ". 


j(X1 +X2) E 
Te 
Il — 
Figure 1. Approximate equivalent circuit of a Figure 2. ZIP model with TAM [22] 


TAM, proposed in the literature [21] 


Where: I; is the current per phase in the stator, Z is the impedance, Rj is the stator resistance, Xm is 
the magnetizing reactance, Io is the current per phase of magnetization, and I is the current per phase in the 
rotor. Aree [23] the equivalent circuit shown in Figure 3 is used to determine, during the execution of the 
load flow, the active and reactive power of the TAM load according to the voltage values calculated in 


Int J Pow Elec & Dri Syst, Vol. 12, No. 4, December 2021 : 2083 — 2094 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 O 2085 


successive iterations. This circuit is transformed by applying Thevenin's theorem [24] and with this, the 

active and reactive power is determined in each iteration. Certain input power is assumed (the work does not 

explain how to obtain this power from the TAM data and its load), and with it and the TAM parameters or 
with added parameters, the slip, and reactive power are determined. Core, mechanical and additional losses 
are not considered. 

Wang et al. [15] a method based on a balance between the active and reactive powers supplied by 
the network and demanded by the TAM is used. From the solution of this balance, the active and reactive 
powers are determined in each iteration of the load flow program. The method assumes that the TAM output 
power is directly proportional to the square of the speed. Arif et al. [3] and Pereira et al. [25] the equivalent 
circuit shown in Figure 3 is also used. This circuit does not consider a parameter to represent the constant and 
additional losses. 

The SimPowerSystem® of MatLab® [26] has a TAM model which requires the parameters and the 

mechanical output power of the TAM. The Power Factory® [27] proposes two methods to represent the 

steady-state model of the TAM in load flow studies: 

- Slip Iteration (Slip Iteration AS). It is based on the exact equivalent circuit. The model equations are 
evaluated in a steady state. The user defines only the EPS input. By performing the iterations 
corresponding to the load flow, the slip and reactive power are determined. 

- Constant PQ load. It is the classical method of representing the TAM. 

Using the equivalent circuit directly as the model of the TAM, as proposed in the existing literature, 

has the following disadvantages: 

- It requires knowledge of all the parameters of the circuit. This information is not provided by almost any 
manufacturer and is complicated when working with several TAMs connected to the same bus, as is the 
case in almost all industries. 

- It does not consider the effect of the type of load or mechanism driven by the TAM. Most TAMs drive 
pumps, compressors, or centrifugal fans, which are variable torque mechanisms whose slip varies with the 
voltage. This aspect is not considered in most of the methods proposed in the literature consulted, nor in 
the SimPowerSystem® or DigSilent PowerFactory® software. 

- They do not consider the effect of frequency variation directly. 

- The models required for power flow studies usually directly reflect the functional relationship of the 
active and reactive power consumed, with the voltage and frequency deviation around its nominal value. 

The consideration of the variation of the frequency and the voltage in the operation of the TAM is an 
aspect that gains relevance in these times due to the increasing use of variable speed drives and of distributed 
generation sources that use non-linear elements [28]. The type of mechanical load that the motors drive is an 
important aspect that must be assessed as it determines their operating characteristics and active and reactive 
power demand from the EPS. About 55% of electric motors drive pumps and fans and 37% drive 
compressors and conveyors [29]. The relationship between speed and mechanical torque of pumps and fans is 
quadratic, and that of compressors and conveyors is constant. In other types of mechanical load driven by the 
motor such as mixers or extruders, the relationship between speed and torque is linear [4]. 


Figure 3. The equivalent circuit is proposed in the bibliography [23] 


In summary, in EPS studies, TAM is represented as a purely electrical load of constant power or 
constant impedance or, at best, expressed by constant slip [17], [30], [31]. In TAM, the behavior of active 
and reactive power consumption to voltage and frequency depends on the type of mechanical load that drives 
and the load factor. In this work, these factors are considered. 


2.2. Proposed model 
Static load models have long been applied to represent static load components, such as resistive and 
lighting loads, and to estimate dynamic load components in EPSs [32]. These models are mainly used in the 
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analysis of the equilibrium condition of EPS [33]. Static load models can be expressed as a polynomial or 
exponential model. The most used is the polynomial model defined as [2]: 


= v y v 
i Po l [p l a t P2 ` Ga) t p| (9 
Pou = Pi * Vou? + Po + Vou + Ps (2) 
v y v 
Q =Q: la: (~~) + q2° (~~) + a| (3) 
Qpu = 91 ` Vou" + q2 ` Vou + q3 (4) 


Where: (pi, p2, p3) are the coefficients of the ZIP model for the active power, Ppu is the active power 
in p.u, Po is the active power at the base voltage, (qı, q2, q3) are the coefficients of the ZIP model for the 
reactive power, V is the voltage, Vbase is the base voltage, Vpu is the voltage in p.u, Qo is the reactive power at 
the base voltage, and Qp is the reactive power in p.u. The parameters (p1, p2, p3) and (q1, q2, q3) must meet 
the conditions: 


Pi +p2+p3 =1 (5) 


qı +42+q3=1 (6) 


As can be seen, these equations have three components: the first, proportional to the square of the 
voltage, corresponds to a constant impedance load, the second, proportional to the voltage, corresponds to a 
constant current load, and the third, as a constant term, corresponds to a constant power load. Because of this, 
this representation of the load is called the ZIP model (Z Impedance, I Current, P Power). 
The equations of the polynomial model to take into account also the frequency variation are [3]: 


P =P |p G +p2-() +p] 04K 4 (7) 
Pou = (Pi + Vou” + Pa Vou + Ps) + (1 + Kpr Af) (8) 
o=o fa GE) +a: () +a] ky af 9) 
Qpu = (41 You” + q2 “Vou +q3) (1 + Kas AP) (10) 


Where: Kyr is the coefficient of variation of the active power, Kg is the coefficient of variation of the 
reactive power, and Af is the deviation of the frequency from the nominal. The model can be represented 
schematically as shown in Figure 4. 


J/g} [1 p2 ps] [q2 q2 q3] 


Figure 4. Schematic of the polynomial model of the TAM 


where f is the frequency of the EPS. 


2.2.1. Estimation of model parameters 
In the application of the polynomial model, an important goal is to estimate the parameters of Figure 
4 (we., pl, p2, p3, q1, q2, q3, Kpf, and Kqf). In this investigation, the exact equivalent circuit shown in Figure 
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5 was used applying the following methodology. With the datasheet or the data provided by the 

manufacturer, the circuit parameters are estimated, following the method explained in [34]. 

- Considering the characteristics of the driven mechanism and using the method explained in the next 
section, the voltage is varied in a suitable range at the nominal frequency. For each voltage value, the 
circuit is solved, and the active and reactive input power is determined. At the same time, the 
characteristic model that relates the powers to the voltage per unit is obtained. 

- The parameters in (2) and (4) are obtained using a curve fitting method. 

- Repeat step 2 but varying the frequency deviation at the nominal voltage. 

- The parameters in (7) and (9) are obtained using a curve fitting method. 


Figure 5. Exact equivalent circuit of the TAM [27] 


where f is the frequency of the EPS. 


2.2.2. Consideration of the influence of mechanical load driven by the TAM 

Next, how the proposed method considers the influence of the type of mechanical load that drives 
the TAM on the behavior of the active and reactive power demand and the voltage and frequency of the EPS 
is explained. The two types of mechanical load that can drive ATMs are: 


— Constant torque loads. 

Tmec = Te (11) 
— Variable torque loads. 

Tmec = To tT en+ T- n? (12) 
where: Tmec is the mechanical torque, Te is the torque at the intersection point of the characteristic curve of 
the motor and the load, (To, Tı, T2) are the variable torque model coefficients, and n is the speed of the rotor 


shaft. The equation of electromagnetic torque produced by the TAM, depending on the parameters of the 
equivalent circuit, is [24]: 


3v1? R 
Tom = ==.: 2 (13) 
7 Os s[(mi+82)' +04 +x2)2| 

where Tem is the electromagnetic torque and œs is the angular synchronic velocity. 
In the steady-state operating zone of the TAM, the following is fulfilled: 

R2 

ee R, (14) 

“2 Sk +X (15) 
Considering (14) and (15), (13) can be written as: 

_ 3M? 
Tem = (16) 


Considering that 3, ws, and R2 are constants, (16) can be written as: 
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The constant K1 can be expressed as: 


ki = Tmecn (18) 


2 
Vin'Sn 


Where K; is the constant used for the electromagnetic torque model as a function of variable voltage and 
constant frequency, Tmecn is the nominal mechanical torque and V ın is the nominal voltage per phase. 
The slip is related to speed through: 


n=n,:(1-s) (19) 


In the case of constant torque load, (11) is used, equating the mechanical torque with that of the load. If the 
voltage decreases, the slip must increase so that the torque remains constant and vice versa. This makes it 
possible to assume that the product (s-V12) remains constant (see (17)) and then: 


k,=V2es (20) 


Where K, is the constant used for the electromagnetic torque model as a function of variable voltage and 
constant frequency. 

The algorithm to be programmed is based on assuming a load torque value for which the 
corresponding Kv is calculated and then varying the voltage from 80% to 120% of the nominal one. For each 
voltage value, calculate the active power and the reactive input power. This makes it possible to obtain the 
characteristics of these two variables as a function of the voltage. 

In the case of variable torque, it should be noted that varying the slip also changes the torque as the 
speed n = ns (1-s) varies. 
where ns is the synchronic speed. 

For the variable torque, it is satisfied that: 


T, =To tT, n+ n? (21) 
Te = To + T,-n,+(1—s) +72 -n,?+ (1-2-5457) (22) 
Te = (To + Ty n + Tz + N57) — (Ty Ng + 2+ heon s+ en. s? (23) 


Substituting (16) in (18) it is obtained that: 


3 


koe (24) 
Equation (17) and (23) and making some transformations are obtained: 

0 = (T+T ns + T2 +n) — (Tins +2- Tn +K vn?) s+ +n? s? (25) 

Ao = (To + Ty -ns + Tz + N57) (26) 

A, = (Tiens +2: Tzn +K- V?) (27) 

A, =T, +n,” (28) 


Where Ao, Ai, A? are variable load and constant frequency model coefficients. 
The equation that gives the slip value for each voltage value is: 
0 =A) — A1 -s HA3: s? (29) 
This equation has two solutions, one of them is within the range of possible slip values (0 <s <1) and 
the other solution exceeds this value by more than 10 times and is therefore not considered. With the slip 
value given by these equations, the equivalent circuit is solved for each of the voltage values and the values 
of active power and reactive power are calculated, as it was done in the constant torque load. 
In the case of frequency variation, in (8) and (10) the voltage is made in per unit equal to 1 and is obtained: 
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Pry = (1+ Kye Af) (30) 


Qpu =(1+ Kar Af) (31) 


Frequency variation implies a new value of synchronous speed and requires a new way of determining torque 
from: 


= 3V1? Ws—Wr 


Tom = oe e (32) 
2 
Tem = qe" (2) - (ws — oy) G3) 


The term V1/œs expressed per unit coincides with the magnetic flux per unit (V/f=k) because it is 
proportional to it. The following term is then defined as “Equivalent Flux”. 


Peq SF (34) 
Substituting (34) in (33) it is obtained that: 
Tem = Ko Peq (ws = wr) (35) 


Where Ko is the constant used for the electromagnetic torque model as a function of variable voltage and 
frequency, Peq is the equivalent flux and œr is the angular shaft speed. The constant Ko is calculated as: 


Ko perns Tmecn (36) 


a 2 
Peqn ‘(@sn-®rn) 


Where eqn is the nominal equivalent flux, @m is the nominal angular shaft speed and @sn is the nominal 
angular synchronic velocity. In the case of constant torque load, frequency values are expressed as: 


f=fa+Af (37) 


Where fn is the nominal frequency. 
This frequency is substituted in: 


w = AEE (38) 


Where p is the number of poles of the motor. 

This equation and (34) with the nominal parameters are replaced in (36). The reactance values are changed 
according to the frequency. With these values and the other TAM parameters, the equivalent circuit is solved, 
and the active power and the reactive power are calculated. If the load is of variable torque, (23) is (35) and 
the speed in rpm (n) is replaced by the speed in rad/s (ws) applying (39). 


new) 
imn Gorr O aje n (8-42 OP 

o?) s+ (r (EY a2) -s? (40) 
Bo = (T +T: (=) -ws + Tp: Sj : ws?) (41) 
B, = (Ka'ba? os +7, (E) s2 n (E) - 08) (42) 
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n- (nC) a 


Where (Bo, Bi, B2) are the variable load and frequency model coefficients. In (40) is then like: 
B —B,:S+B,:s?=0 (44) 


Solving this second-degree equation, two solutions are obtained, the smallest (O <s <1) is valid because it 
corresponds to the normal values of slip in the TAM. The rest of the procedure is the same as for a constant 
torque load. 


3. RESULTS AND DISCUSSIONS 
— Application of the method to a case study 

To illustrate the method, the TAM of the feeding pump of a thermal power plant of 2500 kW; 6000 
V; 2 poles and 60 Hz is selected, whose data is provided by the manufacturer, and parameters calculated from 
these data, appear in [34]. To determine the characteristics of active and reactive power as a function of 
voltage for both constant and variable torque, an operating point at a nominal voltage corresponding to 75% 
of the nominal load was assumed in both cases, at which the active power consumed is P=1817 kW and the 
reactive Q=1134 kVAr. For this, a MATLAB [26] program was used with the algorithm described in a 
previous section. 

Figure 6 (a) shows the characteristic of active power as a function of voltage and in Figure 6 (b) the 
characteristic of reactive power as a function of voltage at nominal frequency. In both cases for a constant 
torque load and a variable torque load. The characteristics vary very little with the type of load. It is only 
noticeable in the active power characteristic, for voltage values are far from the nominal one. 

It is also analyzed how the consumption of active and reactive power varies as a function of the 
voltage when the mechanical load of the TAM also varies, both for constant and variable torque load at 
nominal frequency. Figures 7 (a)-(b) shows the characteristics of active and reactive power as a function of 
voltage. In this case, the torque loads with constant and variable characteristics have the same behavior. 


1.3 1.3 
Constant torque Constant torque 
1.2| — — —Variable torque 1.2 | — — — Variable Torque 
=> Ss 
= d 
a 21.17 
na fot] 
5 = 1} 
=) 
=: a 
20.9 
S 2 
a g 0.8f 
ac 
0.7 
0.6 
0.8 0.9 1 1.1 12 0.8 0.9 1 1.1 1.2 
Voltage (p.u) Voltage (p.u) 
(a) (b) 


Figure 6. These figures are; (a) characteristic of active power as a function of the voltage, (b) characteristic of 
reactive power as a function of the voltage 
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Figure 7. Characteristics of; (a) active power and (b) reactive power as a function of voltage at different load 
factors for constant and variable torque load 
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It is observed that the higher the load, the greater the variation of the active power with the voltage 
and the lower the reactive power. It can also be noted that, in both cases, the lower the load, the linearity of 
the characteristic increases. Using the Basic Fitting option of MATLAB [26], the parameters corresponding 
to the voltage for the model are given by (2) and (4) were determined and appear in Table | for the case of 
variable torque. To determine the characteristics of active and reactive power as a function of the frequency 
deviation for both constant torque and variable torque, an operating point at a nominal voltage corresponding 
to 75% of the nominal load was assumed in both cases, at which the power consumed is P = 1817 kW and Q 
= 1134 kVAr. 


Table 1. Model parameters for variable torque load, at different load factors, and nominal frequency 
Parameters 100% 75% 50% 25% 


pl -0.44 -0.28 -0.13 0.0046 
p2 1.2 074 039 0.13 
p3 0.28 0.53 0.74 0.87 
ql 13 13 12 1.1 
q2 2i A7 1 031 
q3 17 13 078 021 


Figures 8 (a)-(b) shows the active power and reactive power characteristics as a function of the 
frequency deviation in both cases, for a constant torque load and a variable torque load. As can be seen in 
Figure 8 the active power characteristics as a function of the frequency deviation from the nominal, are linear 
with a much greater slope in the case of variable torque loads. However, the variation in reactive power is 
practically negligible. This behavior explains the advantage of using variable speed drives in this type of 
mechanical load driven by the motor for saving electrical energy [8], [35]. Next, the variation of the active 
and reactive power consumption is analyzed as a function of the frequency deviation at nominal voltage, 
when the load factor of the TAM also varies. Figures 9 (a)-(b) shows the variation for constant torque loads 
and Figure 10 (a)-(b) shows for variable torque loads. In both cases, the active power is represented in 
Figures 9 (a), and 10 (a) then the reactive power in Figures 9 (b), and 10 (b). 
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Figure 8. Characteristic of; (a) active power and (b) reactive power as a function of frequency variation 
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Figure 9. Characteristics of; (a) active power and (b) reactive power as a function of frequency variation at 
different load factors for constant torque load 
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Figure 10. Characteristics of; (a) active power and (b) reactive power as a function of frequency variation at 
different load factors for variable torque load 


The active power characteristic as a function of the frequency deviation about the nominal is 
practically independent of the load factor, that is, it can be assumed that it only depends on the type of load. 
Regarding the reactive power characteristic, it can still be assumed that the variation of the reactive power 
with the frequency deviation is practically negligible. Table 2 shows the parameters related to the frequency 
variation of the model of (15). 


Table 2. Model parameters corresponding to frequency deviation, at different load factors, and 
nominal voltage 


Parameters 100% 75% 50% 25% 
Kpf (constant load) 0.81 0.88 0.92 0.93 
Kpf (vatiable load) 2.5 2.6 2.8 2.8 


4. CONCLUSION 

According to the results, the following conclusions can be reached: The representation of the TAM 
through a quadratic polynomial is very suitable for its application in the analysis of EPS performance. The 
most convenient procedure to determine the model parameters is, based on the information provided by the 
manufacturer, to calculate the parameters of the exact equivalent circuit and, through it, to determine the 
active and reactive power characteristics as a function of the voltage and the frequency deviation around 
nominal one. With these characteristics, and applying some curve fitting method, determine the parameters of 
the proposed model. The model parameters depend, in general, on the type (constant torque or variable 
torque) and magnitude of the mechanical load driven by the TAM. The characteristics of active and reactive 
power as a function of the voltage do not depend on the type of load, but on its magnitude, while the 
characteristics as a function of the frequency deviation depend on the type of load but not on its magnitude. 
Active power varies much more with frequency in the case of a TAM that drives a variable torque 
mechanism than if the mechanism is constant torque. The reactive power consumed by the TAM varies very 
little as a function of the frequency deviation. Because of this, it is proposed to disregard it. 
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